Incoherent topological defect recombination dynamics in TbTe 3 
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We study the incoherent recombination of topological defects created during a rapid quench of a 
charge-density-wave system through the electronic ordering transition. Using a specially devised 3- 
pulse femtosecond optical spectroscopy technique we follow the evolution of the order parameter over 
a wide range of timescales. By careful consideration of thermal processes we can clearly identify 
intrinsic topological defect annihilation processes on a timescale ~ 30 ps and find a signature of 
extrinsic defect-dominated relaxation dynamics occurring on longer timescales. 
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Topological defects are non-linear objects which can 
be created any time a symmetry-breaking transition 
occurs. [ll-Q They are predicted theoretically as solutions 
to systems of nonlinear differential equations based on 
Ginzbug-Landau theory. They are of great fundamental 
importance in fields such as cosmology where they ap- 
pear as strings and condensed matter physics where they 
appear in the form of vortices and domain walls. While a 
good understanding of static properties of topological de- 
fects (TD) has come from systems such as liquid crystals, 
the dynamics of TDs are much less understood. Elec- 
tronic phase transitions in charge-density wave systems^ 
are particularly interesting model systems for studying 
the general behavior of the dynamics of topological ex- 
citations. The collective excitations are not overdamped 
which allows the observation of both collective and quasi- 
particle (QP) excitations as they evolve through the tran- 
sition. In particular, they can be used to investigate the 
dynamic behavior of topological excitations such as do- 
main walls in real time using ultrafast laser techniques. 

Recently, time-resolved experiments have shown that 
following a quench created by a strong laser pulse the 
order parameter (OP) oscillates coherently, revealing co- 
herent TD dynamics. [5] Domain walls are created parallel 
to the crystal surface which can coherently annihilate on 
the timescale of a few picoseconds with the accompany- 
ing emission of collective modes which have been detected 
as modulations of reflectivity upon reaching the surface. 
In addition to coherent defect dynamics, incoherent topo- 
logical defects created by the Kibble-Zurek mechanismQ 
are also expected, but very little is known about the dy- 
namics of incoherent TD dynamics in CDW systems, and 
in condensed matter systems in general. 

In this paper we investigate the incoherent evolution 
of TDs using a specially devised 3-pulse femtosecond 
spectroscopy [5] technique which allows the direct back- 
ground free observation of the evolution of the order pa- 
rameter (OP) as a function of time through the elec- 
tronic ordering transition. In a rapid quench experiment 
order emerges in different regions of the sample inde- 
pendently so multiple topological defects can be created. 
Their presence can be detected in the optical response as 



a spatial inhomogeneity of the order parameter. The de- 
termination of incoherent TD dynamics is a challenging 
task, however. Because of thermal diffusion processes, 
which evolve on similar timescales as topological anni- 
hilation and also introduce temperature inhomogeneity, 
careful temperature calibration from independently mea- 
sured frequencies is needed to accurately account for 
thermal effects. We deal with the problem by careful 
calibration of the transient effective temperatures, which 
enables us to unambiguously distinguish the incoherent 
dynamics from thermal diffusion effects. 

In our experiments, we use a three pulse technique de- 
scribed in refs. [f| |g: A "destruction" (D) laser pulse 
at 800 nm excites a cold sample[l3j into the disordered 
state, breaking up the CDW order. We then monitor 
the evolution of the transient reflectivity AR(tp p ) / R ex- 
cited with a weaker pump (P) pulse as a function of 
time delay £dp between the D and P pulse (the pulse 
sequence nomenclature is illustrated in the insert to Fig. 
lb)). The D pulse fluence is adjusted to twice the thresh- 
old for causing the destruction of the ordered state [f|. 
After a quench by the laser pulse, order recovers first 
through the sub-picosecond recovery of the quasiparticle 
gap leading to coherent oscillations of the OP and the 
coherent creation of TDs which decay within 5-8 ps in 
TbTe 3 d. Since the CDW coherence length (-2 nm) 
is much shorter than the size of the laser excited vol- 
ume (~ 50/xm dia), order emerges with different phase in 
different regions, resulting in the formation of topological 
defects whose spatial distribution is determined partly by 
the inhomogeneous excitation and partly by the under- 
lying fluctuations which nucleate the emergence of order 
by the so called Kibble- Zurekf^, [3| mechanism. The re- 
sulting inhomogeneity of the OP leads to observable tem- 
porally resolvable effects in the frequency, linewidth and 
amplitude of the collective amplitude mode, all of which 
are related to the OP, as shown in previous studies|7j. 

In Fig. QJa) we show the raw data on the transient 
reflectivity Ai?(ip p ) / R of TbTe3 at different time delays 
tnp after the D pulse. After the initial QP relaxation we 
observe oscillations due to the coherently excited order 
parameter amplitude mode (AM) and other phonons.0- 
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Figure 1: a) The transient reflectivity of TbTe3 at different 
times £dp after the D pulse. The thin lines are the fits dis- 
cussed in text, b) An example of the FFT power spectra for 
the raw data and the fit at top = 10 ps with and without a 
D pulse. The inset shows the laser pulse sequence. 



IH The level of noise is very small, due to the the excellent 
intrinsic properties of the material, which helps us make a 
detailed quantitative analysis. We analyze the transient 
reflectivity oscillations using the theory for displacive ex- 
citation of coherent phonons:[10] 
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where — (1/r — 7j)/Qj, G(t) = cxp(— 2i 2 /Vp) and Tp 
is the laser pulse length. The first integral represents the 
QP relaxation with the relaxation time r, while the sum 
depicts the response of coherent phonons with frequen- 
cies Qi , and effective dampings jj . A c corresponds to the 
QP relaxation amplitude while the residual value at long 
delays is . To limit the number of fitting parameters 
we keep only two phonon terms corresponding to the AM 
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Figure 2: a) Frequencies of the AM and the 1.7- THz phonon 
as a function of top. b) Effective dampings of the AM and the 
phonon as a function of top- Open squares represent a mea- 
surement from another spot on the sample, c) Amplitudes of 
the AM and the phonon as a function of top ■ The solid lines 
are the frequencies, linewidths and amplitudes calculated us- 
ing an inhomogeneous temperature distribution model(Eq[2]). 



at 2.2 THz and the 1.7-THz phonon which strongly in- 
teracts with the AM at higher temperatures. Q FigQJb) 
shows the fast Fourier transform (FFT) of the raw data 
and of the fit to the data i) without the D pulse and ii) for 
a D-P delay of <dp = 10 ps, clearly showing that Eq. (Q]) 
fits the response very well below ~ 2.4 THz irrespective 
of < D p- 

The ^op-dependence of the frequency, linewidth and 
amplitude of the AM and 1.7 THz phonon modes are 
shown in Fig. The linewidth is shown for two sets of 
data obtained from different spots on the sample, show- 
ing an offset between the two data sets. 

In order to obtain a calibration of the effective temper- 
ature T e ff of the photoexcited sample volume we mea- 
sured independently, by means of a standard pump- 
probe experiment, the T-dependence of the reflectivity 
transients to determine the T-dependent amplitude, fre- 
quency [wam(T)] and damping [7am (J 1 )] for the AM. Us- 
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ing these calibrations we are in a position to determine 
T e g as a function of time from clAm(^dp) and 7am (/dp) 
and take it into account to obtain the thermal inhomo- 
geneity dynamics. The time-dependence of the T e ff is 
shown in Fig. \3j(b). We observe that the two effective 
temperatures T w (iDp) and T 7 (£dp) obtained from the 
wam(^dp) and 7am (tup) systematically differ by approx- 
imately 20 K. This difference arises because the order 
parameter is inhomogeneous within the probe volume, 
which increases the effective linewidth and consequently 
T 7 (irjp)- The order parameter inhomogeneity arises from 
both the thermal inhomogeneity and the presence of topo- 
logical defects. 
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Figure 3: a) The temperature dependence of uiam(T) and 
7am (T) measured in a standard Pp experiment. The in- 
set shows the T— dependence of the amplitudes for the 
AM (circles) and the 1.7-THz phonon (diamonds), b) The 
time-dependence of the effective temperature T e s{tr>p) from 
<^>am(£dp) and 7am(£dp). The solid line is a fit to uam using 
the thermal diffusion model . The inverted triangles corre- 
spond to the calculated T™ M (tDp) taking into account the 
inhomogeneous temperature distribution given by Eq. p]l. 

To account for the thermal inhomogeneity we first fit 
a thermal diffusion model[ll] (TDM) to the effective 
temperature from the AM frequency T u (tnp), and then 
use the TDM parameters to calculate the transient op- 



tical response in the reflectivity which fully takes into 
account inhomogeneity of the temperature in the ex- 
cited volume. As seen in Fig. [3] (b), T w (iop) can be 
fit very well over 5 decades of time from 30 ps to 4 
/is using a one-dimensional ffij TDM, where T(<dp) = 
AT '/ ' \J\ + tDp/ r D + T , and the fit parameter td ~ f20 
ps represents the characteristic heat diffusion timeff f]. 

Next we simulate the optical reflectivity response in the 
presence of the inhomogeneous temperature distribution 
using the thermal diffusion parameters obtained from the 
fit above. The response is given by: 
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where Ai?(ip p )/i?| T ( ZjtDP ) in ((2|) is calculated by interpo- 
lation from the experimentally determined T-dependent 
transients and A op is the optical penetration depth. From 
the response in Eq. ([2]) we calculate the predicted ther- 
mally inhomogeneous linewidth 7am M (^-D-P )> and other 
phonon parameters, as well as the effective temperature 
T' 7 rDM (iDp) in the same way as previously using the tem- 
perature calibration data in Fig. [3] a) . 

Comparing now the simulation with the experiment 
in Fig. [2j we see that the validity of the TDM beyond 
~ 30 ps is well supported by the good agreement of the 
simulation for both the AM and the f .7-THz phonon pa- 
rameters. We can thus safely conclude, that the recovery 
of the order parameter on timescales longer than ~ 30 ps 
is primarily governed by the fD heat diffusion process. 

Below ~ 30 ps however, there is a large discrepancy be- 
tween the calculated T c g , 7am and other phonon param- 
eters in comparison to the data, even after carefully tak- 
ing into account the thermal inhomogeneity. While the 
smaller spatially varying discrepancy for 7am (^dp) [and 
T 7 rDM (^Dp)] at longer £dp might possibly be attributed to 
a spatial variation of the intrinsic AM linewidth, |f 5| the 
observed magnitude and the evolution of 7am(£dp) for 
<dp< 30 ps clearly cannot be assigned solely to the tem- 
perature inhomogeneity. Subtracting the thermal con- 
tribution from the AM linewidth in Fig. 2 b), we can 
now isolate the contribution of the topological defects as 
shown in Fig. @] 

As discussed in the introduction some of the defects 
created in the quench process annihilate coherently re- 
sulting in an aperiodic modulation of the AM intensity 
and frequency in the first ~8 picoseconds The very 
low level of experimental noise in the raw data allows us 
to attribute the observed data scatter to the coherent de- 
fect dynamics in the material, rather than experimental 
noise. 

The modulation of the AM due to the defects that 
annihilate incoherently is unfortunately not detected di- 
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Figure 4: The time-dependence of 7am due to topological 
defect annihilation. Data for spot #1 and spot #2 on the 
sample show significant differences primarily in the long time 
behavior beyond 30 ps, which appears as an offset. 

rectly by our stroboscopic technique. However, the inco- 
herent topological defects give rise to a spatial inhomo- 
geneity of the order parameter and a decoherence of the 
AM oscillations leading to an increased linewidth 7am 
for <dp < 30 ps which we have detected in our experi- 
ments. Concurrently the defects give rise to a softening of 
the collective mode wam, because of the OP suppression 
which they cause. The increase of the coupled 1.7-THz 
phonon effective damping at shorter <dp as shown in Fig. 
[5] is presumably also caused by the inhomogeneity of the 
OP. A further manifestation of the incoherent annihila- 
tion is the increase of the amplitudes of the AM and the 
phonon with respect to the TDM for iDP <^ 30 ps, which 
is also consistent with the suppression of the OP when 



one takes into account that the amplitudes of the modes 
increase with the decreasing OP amplitude deduced from 
their T-dependence shown in the inset of Fig. [3ta). 

Apart from intrinsic topological defect annihilation 
processes, which we have identified on a timescale of 30 
ps, we expect to observe annihilation of domain walls 
pinned to defects and imperfections at longer times. The 
timescale of their annihilation may extend well beyond 
30 ps into microseconds. Evidence for such slower ex- 
trinsic recombination processes conies from the spot-to- 
spot variation shown in Figs. and 4. The dynam- 
ics on the shorter timescales up to 30 ps appears qual- 
itatively the same in different spots, but the dynamics 
at longer timescales appears as a background which is 
strongly spot-dependent, as expected for pinned defect 
annihilation dominated by extrinsic sample imperfections 
or impurities. 

In conclusion, these experiments demonstrate the pos- 
sibility of studying both coherent and incoherent topolog- 
ical defects dynamics in complex materials in which the 
order parameter can be monitored in real time through 
the dynamics of the collective mode. The dynamics on 
a timescale of ~ 30 ps can unambiguously be associated 
with intrinsic topological defects annihilation in TbTe3 
following a laser quench arising from the time-dependent 
inhomogeneity and suppression of the order parameter. 
The inhomogeneity causes an increased effective damp- 
ing of the amplitude mode while the suppression of the 
order-parameter is indicated by an additional softening of 
the AM-mode frequency. Beyond ~ 30 ps we find a pre- 
dominantly thermal-diffusion governed order-parameter 
dynamics with a signature of extrinsic defect annihila- 
tion dynamics, suggested by significant variation from 
spot to spot across the sample. 
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The growth of the samples by a self-flux technique and 
subsequent characterization was described elsewhere. [L?T] 
The diameters of the beams are much larger that the 
optical penetration depth so on the relevant timescale 
the heat diffusion is ID. 

In the standard 2-pulse pump-probe experiments we ob- 
serve some spatial variation of the AM linewidth which 
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we assign to the intrinsic structural-defect/impurity den- 
sity spatial variation in the sample. Although the care 
was taken that T-dependent standard 2-pulse Pp scans, 
used to obtain the data for simulation, were performed 
on the same spot on the sample (spot #1) as the 3-pulse 
DPp scans, some small daily drift of the spot position 
on the sample can not be excluded. A possible contri- 



bution to the AM linewidth discrepancy at long tup by 
the spatial variation of the intrinsic AM linewidth is sug- 
gested by the larger discrepancy, observed beyond ~ 30 
ps, between the data from another spot (spot #2) and the 
corresponding TDM simulation based on the first spot 
(spot#l) T-dependent scans [see Fig. [2](b)]. 



